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The aim of this note is to bring attention to a simple class of discrete dynamical 
systems exhibiting some complex behaviour. Each of these systems is defined as 
a self-mapping of the unit square I 2 = I x I (with / the closed unit interval [0, 1]) 
and is obtained by coupling two families of self-mappings of the interval I. There 
is no real mathematics to be found here (in the sense of results stated and proved) 
and in fact there is an almost complete lack of precise statements. The only thing 
on offer is the definition of the mappings and a few nice pictures showing examples 
of their asymptotic limit sets. I have written a JavaScript program, accessible at 
www.math.uni-bielefeld.de/~preston/iterates.html, which can be used to 
'discover' more about these mappings. The program might prove to be helpful 
for anyone interested in doing this. 

A mapping h : Z — > Z of a set Z into itself will be considered as a discrete 
dynamical system in the usual way: The set Z is the state space and if the 
system is in state z at time n then h(z) is the state at time n + 1. For each z G Z 
the sequence {z n } n > Q , where z — z and z n+x = h(z n ) for all n > 0, then describes 
the successive states of the system given that the system started in state z at 
time 0. This sequence is called the orbit of z under h. 

For various classes of mappings it is expected (or hoped) that some kind of limit 
set should exist. Assuming Z is a metric space with metric d this means that 
a closed subset S of Z should exist such that the orbit {z n } n >o of each 'typical' 
point z converges to the whole of S. An orbit {z n } n >o converging to the whole 
of S means that the following two statements hold: 

- The orbit eventually comes arbitrarily close to the set S. More precisely, for 
each e > there exists m > such that d(z n , S) < e for all n > m (i.e., for 
each n > m there exists y n G S with d(z n , y n ) < e). 

- Each y G S is an accumulation point of the orbit. More precisely, for each 
e > there exist infinitely many indexes n with d(y, z n ) < e. 

(Note that if there is a closed set S satisfying these two conditions then it is 
uniquely determined by the orbit {z n } n >Q.) That this should hold for a 'typical' 
point z means that the set of points for which it does not hold should be negligible 
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in either a measure-theoretical or topological sense (i.e., it should have measure 
zero with respect to an appropriate measure on Z or be of the first category in 
the terminology of the Baire category theorem). 

Suppose that a limit set S is thought to exist for some explicitly given mapping 
h : Z — > Z. Then an approximation to S can be displayed on a monitor as follows: 
Start with some random initial point z (which it is hoped will be 'typical') and 
consider the orbit {z n } n >o of z under h. Compute the first N terms of this orbit, 
where N is large enough so that z n is less than one pixel from S for all n > N and 
then display the next M terms of the orbit on the monitor, where M is chosen so 
that the pixels of the displayed points more-or-less fill out the set S. Of course, 
these statements concerning the choice of M and N are extremely vague and in 
practice it is necessary to resort to trial and error to find suitable values for them. 

Conversely, suppose it is not known whether a limit set exists, but when the 
above procedure is carried out for many different initial points z and for ever larger 
values for N the resulting image is always the same. Then this can be interpreted 
as evidence that a limit set exists and that the image is an approximation to the 
limit set. In this sense there is strong evidence that each of the two-dimensional 
dynamical systems which we now introduce possesses a limit set. 

As already stated, the dynamical systems we are interested in are defined in the 
unit square I 2 and are obtained by coupling two families of self-mappings of the 
interval I. Each such family of self-mappings of I will be given by a mapping 
/ : I 2 — > I, where for each p G I the mapping x h- > f(p,x) is the self-mapping 
of / corresponding to the parameter value p. We write f p (x) instead of f(p,x) 
to emphasise that the first argument is to be regarded as a parameter. The 
prototypical example is the logistic map (or family) £ with 

£ p (x) = Apx(l - x) . 

Another typical example is the tent map (or family) t with 

t p (x)=p(l - \2x-l\) . 

As well as two families of self-mappings of / we also require two couplers. A 
coupler is just a mapping c :/—)■/, for example there a simple linear coupler 
with c(x) = b + rx, where b, r are constants with b > 0, r > and b + r < 1. 
This coupler will be referred to as the linear+ coupler and the parameters b and 
r as its base and rate. 

Given two families of self-mappings f,g:I 2 ^>-I and two couplers c, d : I —> I 
the coupled system is then the mapping h : I 2 — )■ I 2 defined by 



h(x,y) = (f c{y) (x),g d{x) (y)) 
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for all (x,y) G I 2 . The first component of h(x,y) is thus obtained by applying 
the family / to the first argument x with the parameter value given by applying 
the coupler c to the second argument y. In the same way the second component 
of h(x, y) is obtained by applying the family g to the second argument y with the 
parameter value given by applying the coupler d to the first argument x. 

To be more explicit, in the following let us take both / and g to be the logistic 
map i and c and d to be linear+ couplers, c with base b and rate r and d with 
base b' and rate r'. Then the mapping h is given by 

h(x, y) = (4(6 + ry)x(l — x),4(b' + r'x)y(l — y)) 

and so h depends on the four parameters b, r, b' and r'. Carrying out the 
procedure described above for this mapping h with b — b' — 0.4 and r = r' = 0.6 
results in the following image: 




The values iV = 1000000 and M = 100000 were used here and the initial value 
was (0.7, 0.6). Changing the initial value or increasing iV does not alter the image. 
The value for M seems about right, but this is rather subjective and also depends 
on the monitor. 

Now what is more interesting than the limit set for a particular mapping is how 
the limit set behaves as a function of the parameters b, r, b' and r'. If we choose 
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a line segment (or more general curve) in this four- dimensional parameter space 
and compute the limit sets along a grid of points on the curve then we obtain a 
succession of images which can be made to appear as a 'video'. 

Consider the simple case in which 6' and r' are fixed with b' = 0.4 and r' = 0.6, 
where r = 1 — b and b runs from to 1. In other words we are considering the 
fine segment si—)- (s, 1 — s, 0.4, 0.6) in parameter space. Here is a selection of 
limit sets from points lying on this line segment: 




6 = 0.00 r = 1.00 6 = 0.05 r = 0.95 6 = 0.10 r = 0.90 




6 = 0.15 r = 0.85 6 = 0.20 r = 0.80 6 = 0.25 r = 0.75 




6 = 0.30 r = 070 6 = 0.35 r = 0.65 6 = 0.40 r = 0.60 
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6 = 0.60 r = 1.85 




6 = 0.75 r = 0.25 6 = 0.80 r = 0.20 6 = 0.85 r = 0.15 




6 = 0.90 r = 0.10 6 = 0.95 r = 0.05 6 = 1.00 r = 0.00 



The images shown above are taken from a uniform grid of points along the line 
segment and give a very rough impression of some of the different kinds of limit 
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sets which can occur. (The images where the limit set is a finite periodic cycle 
have been manipulated by making the points larger. Before this was done they 
were rather hard to see.) In order to get a more complete picture, however, it 
is necessary to consider a much finer grid of say 1000 points (which can be done 
using the JavaScript program). 

The qualitative behaviour of the limit sets along any line segment of the form 
s (->■ (s, 1 — s, b', r') with b' + r' — 1 seems to be essentially the same as that along 
the line segment s (->■ (s, 1 — s, 0.4, 0.6). Moreover, if one or both of the logistic 
maps is replaced with the tent map and the limit sets are considered along one 
of these line segments then the resulting succession of images is basically similar 
to that from the original example. 

We now look at something slightly different. Again consider two families of self- 
mappings f,g:I 2 -^-I and two couplers c, d : I —> I. Then there is another way 
of obtaining a coupled system: This uses the mapping hi : I 2 — >■ I 2 defined by 

h'(x,y) = (fc(y)(x),gd(f c(y) (x))(y)) 

instead of using h. To see what is going on here let us write the definitions of h 
and hi one under the other in the following form: 

h(x,y) = (x',y), where x = f c ( y )(x) and y = g d(x )(y), 

h'(x,y) = (x',y'), where x' = f c ( y )(x) and y' = g d { x '){y)- 

From this it can perhaps be seen how the mapping h' could be 'discovered' as a 
result of making a common programming mistake when writing a program for h. 

The reason for looking at this new class of mappings is that there is again strong 
evidence that the mappings possess limit sets, but the behaviour of these limit 
sets along the usual line segments seems to be very different to that which was 
seen above. As before, let us restrict our attention to the case in which both / 
and g are the logistic map t and c and d are linear-l- couplers, c with base b and 
rate r and d with base b' and rate r'. Then the mapping h! is given by 

h'(x, y) = (4(6 + ry)x(l - x),4(b' + 4r'(b + ry)x(\ - x))y(\ - y)) . 

Consider the same line segment s h- >■ (s, 1— s, 0.4, 0.6) in parameter space as above, 
i.e., with V and r' held fixed and taking the values 0.4 and 0.6 respectively, with 
r = 1 — b and with b running from to 1. Here are the limit sets taken at the 
same uniform grid of points along the line segment: 
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b = 0.45 r = 0.55 b = 0.50 r = 0.50 b = 0.55 r = 0.45 
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6 = 0.60 r = 1.85 6 = 0.65 r = 0.80 6 = 0.70 r = 0.75 




6 = 0.75 r = 0.25 6 = 0.80 r = 0.20 6 = 0.85 r = 0.15 




6 = 0.90 r = 0.10 6 = 0.95 r = 0.05 6 = 1.00 r = 0.00 



(Again, where the limit set is a finite periodic cycle the points have been enlarged 
to make them easier to see.) As in the first case, the qualitative behaviour of the 
limit sets along any of the lines segments considered here seems to be essentially 
the same as that along the line segment s i— > (s, 1— s, 0.4, 0.6). Moreover, replacing 
one or both of the logistic maps with the tent map does not change the resulting 
succession of images beyond recognition. 
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